Unconventional pairing in the iron pnictides 
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We determine the anisotropy of the spin fluctuation induced pairing gap on the Fermi surface 
of the FeAs based superconductors as function of the exchange and Hund's coupling Jh- We find 
that for sufficiently large Jh, nearly commensurate magnetic fiuctuations yield a fully gapped s - 
pairing state with small anisotropy of the gap amplitude on each Fermi surface sheet, but significant 
variations of the gap amplitude for different sheets of the Fermi surface. In particular, we obtain the 
large variation of the gap amplitude on different Fermi surface sheets, as seen in ARPES experiments. 
For smaller values of Hund's coupling incommensurate magnetic fiuctuations yield an s -pairing 
state with line nodes. Such a state is also possible once the anisotropy of the material is reduced 
and three dimensional effects come into play. 
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I. INTRODUCTION 

The recently discovered FeAs based family^ has been 
captivating the community primarily because of its high 
superconducting transition temperatures, with Tc values 
vifell above 50K in some casesi^i^i^ While such values for 
Tc could potentially be due to the interaction between 
electrons and lattice vibrations, the vibrational modes 
of the common structural unit, the FeAs -planes, are 
rather low, making electron-phonon interactions as the 
sole or primary mechanism unlikely.— The observation of 
antiferromagnetic order in undoped systems at ambient 
pressure^ has therefore been one of the key motivations 
to explore spin fluctuations as the primary mechanism 
for superconductivity in the pnictides.— i^°i^^ In this case, 
the role of phonons, as intermediate boson and pairing 
glue, is being played by collective paramagnon excita- 
tions of the electron fluid. In order to determine which 
many body interaction is responsible for the formation 
of Cooper pairs, an understanding of the symmetry and 
detailed momentum dependence of the pairing gap is cru- 
cial. 

Experimentally, the strongest indication that the pair- 
ing gap in the pnictides has line nodes comes from nuclear 
magnetic resonance (NMR) measurement with power 
law variation of the spin lattice relaxation rate, T^^ oc 
2^3 12,13,14,15 Qj^ ^jjg other hand angular resolved photoe- 
mission spectroscopy (ARPES) experiments flnd node- 
less, weakly anisotropic gaps on the Fermi surface ji^^ii^ 
Penetration depth measurements in the 122-compound 
Bao.93Coo.o7Fe2As2 support gap nodes,~ while mea- 
surements for the 1111 system NdFeAsOo.gFo.i favor 
anisotropic gaps that remain flnite everywhere on the 
Fermi surface.— Interestingly, NMR results of Ref.— and 
the ARPES data of Refs.— ^^^ are consistent to the extent 
that they see evidence for multiple gap values. ARPES 
measurements demonstrate that the two Fermi surface 
sheets around the F point of Bao.6Ko.4Fe2As2 have am- 
plitudes that differ by more than a factor 2M' Knight 
shift and spin lattice relaxation rate measurements in 
PrFeAsOo.sgFo.iiwere fit to two gaps with ratio ~ 3.2;^^ 

In this paper we determine the momentum depen- 



dence of the superconducting gap, where Cooper pair- 
ing is due to the exchange of antiferromagnetic spin fluc- 
tuations. We flnd, in agreement with previous calcula- 
tions, Ref.— >i2ii2ii2^ that the pairing symmetry is extended 
s-wave with the gap on different Fermi surface sheets 
being out of phase, i.e. we flnd an s^ pairing state. 
Superconductivity is caused by the enhanced collective 
spin-fluctuations in the proximity to an ordered anti- 
ferromagnetic state. We flnd that commensurate mag- 
netic correlations can be caused by including a sufflciently 
large Hund's rule coupling Jh, even in an itinerant mag- 
netic material. We also flnd that a large Hund's cou- 
pling generally yields a stronger tendency towards su- 
perconductivity where transition temperatures of 50K 
are possible. We demonstrate that the gap function is 
weakly anisotropic for most sheets of the Fermi surface, 
while a signiflcant anisotropy remains. Depending on the 
strength of the exchange and Hund's couplings Jh the 
gap of this Fermi surface sheet vanishes on Hne nodes 
(for small Jh) or exhibits a moderately anisotropic vari- 
ation along the Fermi surface (for larger, more realistic 
values of Jh)- We also comment on the fact that a sizable 
interlayer coupling, as relevant for the 122 FeAs-family, 
might lead to a nodal superconducting state while for 
the more anisotropic 1111 family a fully gapped state is 
more likely. A possible explanation for the conflicting 
ARPES and NMR flndings is that experiments sensitive 
to the maximum of the gap, such as ARPES, see large 
gaps, while experiments sensitive to the minimum of the 
gap, such as NMR, flnd node Hke features due to im- 
purity induced states in the gap.-S^ The latter is due to 
the fact that non-magnetic impurities in an s^ pairing 
state behave like pair breaking magnetic impurities in a 
conventional s-wave superconductor. 

The spin fluctuation approach relies on two key 
assumptions:^ i) the proximity to a magnetic instabil- 
ity with paramagnons as relevant collective modes, and 
ii) conventional Fermi liquid behavior away from the in- 
stability. While electronic correlations of the Fe 3-d or- 
bitals in the pnictides are relevant, the multi-orbital na- 
ture of the system is Hkely the reason that strong local 
correlations, reminiscent of a system close to a Mott in- 



sulating state do not seem to be dominating. In addi- 
tion, the carrier density of the FeAs systems does not 
seem to be anywhere close to an odd number of electrons 
per Fe 3-d site, strongly suggesting that there are no 
Mott-Hubbard bands with appreciable spectral weight. 
Rather, these systems are closer in their behavior to 
band-insulators or semimetals, however with the bot- 
tom of the electron band somewhat below the top of a 
hole band. The latter leads to the observed hole and 
electron sheets of the Fermi surface. Consistent with 
this picture is that undoped ambient pressure systems ex- 
hibit a small but well established Drude conductivitj*^ 
and magneto-oscillations^ in what seems to be a par- 
tially gapped metallic antiferromagnetic state. Above the 
magnetic ordering temperature a sizable Drude weight, 
not untypical for an almost semimetal has been observed. 
The magnetic susceptibility of BaFe2 As2 single crystals^^ 
above the magnetic transition is only very weakly tem- 
perature dependent and shows no sign for local moment 
behavior of the Fe 3-d electron spins. AT-ray absorption 
spectroscopy for LaFeAsOi_a;Fj; is consistent with a rigid 
band filling upon F-doping and moderate values for the 
effective Hubbard interactioni^ Clearly, these observa- 
tions do not imply that the interactions in the FeAs sys- 
tems are weak, but rather that the phase space for strong 
local correlations is limited and suggest that predomi- 
nant electron-electron interactions are related to inter- 
band scattering between the hole and electron sheets of 
the Fermi surface. Despite very interesting approaches 
based upon the assumption that the FeAs system are 
doped Mott insulators ;22422 we take the view that the 
iron pnictides may be good examples for a system where 
collective longer ranged spin and charge excitations play 
an important role. 

As shown first by Berk and Schrieffer^ magnetic fiuc- 
tuations suppress pairing for a gap function ^aia2 (p) = 
Aq that is constant as function of momentum p and band 
indices a^. However, changing the sign of Aa^aa (p) as 
function of either p or ai , 02 allows for nontrivial super- 
conducting states due to paramagnon fluctuations and 
makes such fluctuations a powerful pairing mechanism. 
In case where only one band contributes to the Fermi 
surface the sign change is a function of momentum p, 
and may lead to line or point nodes of the gap. If there 
are several bands crossing the Fermi energy, strong in- 
terband scattering can lead to a sign change of the gap 
between different Fermi surface sheets, without leading 
to gap nodes. The s^ -pairing state that results from 
our analysis was proposed in the context of the FeAs 
systems in Ref.— in a model with structureless (in mo- 
mentum state) interband pairing interactions. In such a 
state, one would always obtain fully gapped Fermi surface 
sheets. Our analysis shows that the model of Ref."? cap- 
tures the s^ state properly but that one needs to include 
the momentum dependence of the pairing interaction to 
obtain states with residual anisotropy of the pairing gap, 
including states that possess nodes of the gap on a given 
Fermi surface sheet. A careful investigation of the role of 



interband scattering in systems with close to perfect nest- 
ing between distinct Fermi surface sheets was performed 
in Refs.— and^^. These approaches demonstrate that 
under certain circumstances, pairing interactions are en- 
hanced due to interband nesting. At the level of the weak 
coupling expansion used in Ref.— , this conclusion does 
depend on whether the pairing mechanism is due to spin- 
orbital or charge fluctuations. Our results are consistent 
with these flndings, but favor a spin-fluctuation mech- 
anism boosted by intrasite, and inter-orbital exchange 
and Hund's rule coupling. Our approach is closest to the 
results of Refsi^ii^. The key emphasis in our work, as 
compared to these interesting investigations, is to quan- 
titatively analyze the variation of the pairing gap on in- 
dividual Fermi surface sheets as well as between distinct 
sheets. 



II. THE MODEL 

Electronic structure calculations clearly show that the 
states close to the Fermi level are predominantly of Fe- 
3d character with several sheets of the Fermi surface,— 
as conflrmed in recent angular resolved photoemission 
spectroscopy (ARPES) experiments i^^i^^i^^i^^ Given the 
need to change the sign of the gap function l^aia-i (p), 
this leads to the proposal by Mazin et al^ that the gap 
function on sheets coupled by the magnetic wave vector 
are out of phase. 

We use a tight binding description of the Fe-da;^, dyz 
states of the FeAs systems identical to the one proposed 
by Raghu et alM- There are two Fe atoms per crystal- 
lographic unit cell leading to the tight binding Hamilto- 
nian: 



H^ 



E 

p,a/3,cr 






p/3<T 



(1) 



pQ,^ is the creation operator of an electron with 



where dj, 

momentum p and spin a. a refers to the orbital degree 
(i.e. xz and yz) as well as the label of the Fe atom within 
the unit cell. Momenta go from — 7r/a to it /a where a — 
V2ao with Fe — Fe distance hq. Thus Ep is a (4 x 4) 
matrix. As in Ref.— we assume, for simplicity, that all As 
atoms in the unit cell are identical. This approximation 
seems justifled as there are virtually no As states close to 
the Fermi level. The primary relevance of the As states is 
only to determine the indirect overlap between Fe orbitals 
on different sites. With these assumptions, we obtain a 
block structure for the tight binding Hamiltonian of the 
form 



E-p = hp i 



1 



(2) 



with (2 X 2) unit matrix 1 and Pauli matrix t^. hp is a 
diagonal (2 x 2) matrix with diagonal elements 

h^ — 2^2 cos {pxa) + 2^3 cos {pyu) 

2t3 cos (pxtt) + 2t2 cos (pya) (3) 



hi' 




the spin density wave Q = (f : ■^) becomes Q = ( 0, ^ ) 
in the larger BZ. 

In Fig. [T]we show the Fermi surface that results from 
the above tight binding parametrization at a density n — 
1.05. To illustrate the two Brillouin zones used in the 
above discussion we plot the Fermi surface in an extended 
zone scheme. To make contact with Ref.— , we note that 
the axis defining the dxz and dyz orbitals are rotated by 
7r/4 relative to each other. 

Next we include the local electron-electron interaction 
into our theory and write 



Hint = U^nia^fliai+U' ^ ?^^a'^^6 
2, a i,a>b 

-Jh 2_j ( ^®*a • Sib + - 



i,a>b 



Z^ia^^ib 



iba 



(6) 



Figure 1: (Color online) Fermi surface of the tight binding 
parametrization described in the text in the Brillouin zone 
that corresponds to two Fe atoms per unit cell (red diamond 
with axes labeled by k^ and ky) and the larger Brillouin zone 
that corresponds to a unit cell with one atom per unit cell 
and the larger Brillouin zone that corresponds to a unit cell 
with one atom per unit cell (solid square with axes labeled by 
Px andpy), respectively. 



Both diagonal elements of the (2 x 2) matrix 5^ are 

6l^-Sf=At,co.{^)cos{^^) (4) 

while the off-diagonal elements are 



rl2 r21 a, ■ (VxO' 

'p ='^p =4f6sm(^ — 



sm 



PyO. 

2 



^ (5) 



where riiacr = d\^„diacr is the occupation of the orbital 
a with spin a at site i. Uia = J2a "^io-" ^^ ^^^ total 
charge in this orbital and s^a = \ Y^aa' d\aa'^<y<y'diaa' the 
corresponding spin. Thus, we include intra- and inter- 
orbital direct Coulomb interactions, U and U' as well as 
the Hund's rule coupling Jh and the exchange interac- 
tion J. The latter are of interest as they affect the spin 
correlations of electrons in different orbitals. In what 
follows we use U — leV, U' = 0.5eV, electron den- 
sity p = 1.05 per site, and we vary J = Jh between 
J = and J — 0.5eV to explore the role of the exchange 
and Hund's interactions on the pairing state. Recent X- 
ray absorption spectroscopy measurements support val- 
ues for the Hund's couplings that lead to a preferred high 
spin configuration,— leading to larger values of Jh ^ U. 
The interaction term can be put into a more compact 
form^^ 



The individual parameters, determined from fits to full 
potential density functional calculations for LaFeAsO are 
t2 = 0.495eV, h = -0.026eV, h = -0.026eV, h = 
-0.36 eV. 

Because of the assumption of treating all As atoms 
identically, regardless of whether they are located above 
or below the Fe planes, we can describe the system in 
a unit cell with only one Fe atom and can upfold the 
band structure into a larger Brillouin zone, i.e. we ob- 
tain a (2 X 2) matrix tight binding £k in the larger Bril- 



louin zone. It holds hp ~ h 
with reciprocal lattice vector G 



a eic 



; and (5p+G = -5^ 
^ , O) and we obtain 

i5p for states in the original, smaller Bril- 



(5p for momenta outside of 



louin zone and ek+cj = h-p 

it. The momentum k in the new, larger Brillouin zone, 

with -^ < kx.v < :?-, is given hy k^ = ^ {px - Py) 



a.0 



ag ' 



and ky ~ -^ {px +Py)- For example the wave vector of 



Vj, 



V2 



Hint — T 



rjaiaa.agai jt jt J. J. (y-) 



and, in the absence of spin orbit interaction, split into a 
spin and a charge contribution: 



T-7-aia2,a3a4 

^0-i(T2,0-3<''4 



:C/° 



^aia4 ' ^(72 era 



n c "cricr4"o-2 



(8) 



The above Hamiltonian is then recovered if we chose 

U if Oi = 02 = 03 = 04 
Tjaia4,a2a3 ^ ) U' if Oi = O3 7^ 02 = O4 , , 

Jh if oi = 04 7^ 02 = 03 ^ ' 
J if Oi = 02 7^ 03 = 04 



for the spin part of the interaction, and 



with 



U" 



U if ai = 02 = 03 = 04 

— [/' + 2Jh if ai = 03 7^ 02 = 04 

2U' — Jh if fli = 04 7^ ^2 = 03 

J if ai = 02 7^ (23 = 04 



(10) 
for the corresponding charge contribution, respectively. 



A. Collective spin and charge fluctuations 

We determine the single particle and collective mag- 
netic excitation spectrum within a self consistent one 
loop approach, the multiple orbital version^^*^ of the 
fluctuation exchange approximation of Ref.— . Once we 
have self consistently determined the fermionic Green's 
function G"'' (fc) where k = (k, Wn) stands jointly for 
the crystal momentum k and the Matsubara frequency 
ujn — (2n + 1) ttT, we determine the symmetry of the 
pairing state from the linearized gap equation. In the 
normal state, the matrix Green's function of the prob- 
lem is 



G (k) = {iujj - Ek - S (fc) 



(11) 



where Gk, Sfc, £k are all 2 x 2 matrices in orbital space 
in the larger Brillouin zone. The self energy is given as 
a sum of a Hartree-Fock contribution and a fluctuation 
term 






W = 1.1.G' 



.a3a4(;./)pa^a3.a4a2^^_^/^ (12) 



where '^k ■ ■ ■ ~ W^ Sk n • ■ • includes the summation over 
momenta and over Matsubara frequencies. 

Introducing the particle quantum numbers A = 
(01,02) and B = (03,04) labeling the rows and columns 
of two particle states interaction, y'^^"'^'""^"'^ (q) = rf^ (q) 

becomes a 4 x 4-dimensional symmetric operator Fp/j (q). 
Similarly we obtain in this two particle basis a matrix 
representation for the spin and charge couplings C/* and 
U": 



U' 



and 



U" 



U Jh 

U' J 

J U' 

Jh U 



U W 

W J' 

J' W 

W U 



(13) 



(14) 



where W = 2U' - Jh and W = 2Jh - U' . In this 
two particle formalism it is now straightforward to sum 
particle-hole ladder and bubble diagrams and it follows 



1- 



^Ph (q) = t:V' (q) + t:V' (?) 



(15) 



V^iq) = U^[l-x{q)U^] x{qW (16) 

1 



:U^x{qW 



v 



(q) = U-{l + x{q)U^) 'x{q)U' 



-U''x{q)U' 



(17) 



Here x (?) is the matrix of particle-hole bubble in the two 
particle basis. Explicitly it holds: 



T 



X 



aia2,a,a^ ^^-^ ^ _ J^G"^"'' {k + q) G""""' {q) . (18) 



The Hartree-Fock term of the self energy 






\ '' I _TJs,o.3ai,aia2 7^0,0301,0402 j 7^0304 /_-\ 

0304 ^ ' 

(19) 
is frequency and momentum independent and determined 

by Gq'"''' (r") = (4a3'^oo4). It holds for the diagonal 

elements 

S^T = Una, + {2U' -Jh)J2 "'^^ 

02^01 

= {U-2U' + JH)na,+{2U'-JH)n (20) 
whereas the off-diagonal elements (oi y^ 02) are given as: 

S^^;^ = (2 Jh + J - U') (dl,^ do,a. ). (21) 

We are interested in the superconducting transition 
temperature and the symmetry of the superconducting 
state, determined by the corresponding anomalous self 
energy $k {'^n)- Summing up the same class of diagrams 
in the superconducting state yields 

$°i°2(fc)= ^ Tll''^'''-"'^k ~ k') F"--""^ {k') , (22) 

/i;'a304 

with Gor'kov function F (fc). Fpp {q) is the corresponding 
operator in the two particle representation. In this paper 
we only solve the linearized version of Eq. l(22|) to deter- 
mine the superconducting transition temperature as well 
as the nature of the pairing state right below Tc. Close to 
the superconducting transition temperature we linearize 
the anomalous propagator 



F{k)~~G{k)^{k)G{-k) 



(23) 



and obtain 
$0102 (-^) 



T 



Y^ p030i,0204 (^ _ ;,/) 

k' a-^a^^a^aQ 
Qa^a, (^f,f^ ^0506 (^/) ^0604 (•_^/) _ (24) 



Since F"^"* (fc') is of first order in the anomalous self 
energy $'*i'*2 (^k), the Hnearized version of Eq. (|24l) is 
determined by f pp (g) for $°i"2 (fc) ^ q. In this limit 
it follows, after summing the same bubble and ladder 
diagrams as for Tph (q) that 



with 



1; 



^P (.1) = 2^' ('^) - 2^' ^'^^ 



(25) 



the group Dih, i-e. 

£k = ^£k for all R e Dill 



(27) 



Since the two orbitals d^z and dyz transform like coordi- 
nates for in-plane symmetry operations, it holds 



Rel, = D 



(1)? (n^-^A ^ 



(28) 



V^{q) = U^{l~x{q)U') 'x{q)U^ + ^ 

V^q) = U-{l + xiq)U^y'x{q)U^ -^- (26) 

In what follows we first solve the coupled equations 
Eqs. (HH), Jul), (dSl) and JUl) in the normal state on 
a 32 X 32 lattice with 2^^ Matsubara frequencies. The 
solutions of the normal state equations are then used to 
solve the linearized equation for the superconducting self 
energy. In order to determine the superconducting transi- 
tion temperature we replace $"i°2 ^^^ q^ ^j^g j ^ g of Eq. 
(|24|) by A^*^^"^ (p). The resulting eigenvalue equation 
yields an eigenvalue A = 1 if T = Tc, i.e. the temperature 
where the linearization is permitted. For T > Tc, it holds 
A < 1 for the largest eigenvalue. Even if A < 1, the result 
is still useful as (1 — A)~ is proportional to the pairing 
correlation function. Most importantly, the eigenvector 
of the leading eigenvalue determines the momentum and 
band-index dependence of the gap right below Tc. In 
order to simplify the above eigenvalue equation we re- 
place Fpp (p) by its zero Matsubara frequency value, i.e., 
Fpp (p, w„ — 0). Thus, we keep the dynamic excitations 
that determine the frequency dependence of the normal 
state single particle self energy, but assume that the dy- 
namics of the pairing interaction is structureless. Such an 
approximation would be problematic close to a magnetic 
quantum critical point with diverging antiferromagnetic 
correlation length;^ but is expected to be reasonable for 
intermediate magnetic correlations, as seems to be the 
case in the FeAs systems. A consequence of this ap- 
proximation is that we lose the information about the 
frequency dependence of the anomalous self energy. We 
keep its momentum and orbital index dependence. 



where D^^ is the representation of R which transforms 
the coordinates. It then follows that the spinor 



Cka 



Ck,yz,Cr 



transforms as 



Rck^D'}^^ 'c^(i,^ 



(29) 



(30) 



which determines the transformation properties of the 
superconducting gap function in the singlet channel: 

i?$f = Y. D^^l;,'D^^l^'^if,,^ . (31) 



a'h' 



It follows for the transformation of the gap under the 
point group operations: 



E^ 



yrZx ,Ky ) 



^ ^."'a: T"'-y / 



'^ (."•a: ■t'^y ) 









ni 



yKy J Kx) 

[_ f^x 7~~ "^y ) 
\ i^x 1 "'y / 



\kj^ . — ky) 



ky) 



^yy ^yx 

^xy (^xx 

i^Ky ,Kx ) yi^y ^i^x ) 




(32) 



B. Symmetry considerations 

For a proper interpretation of the momentum depen- 
dence of the superconducting gap in a multi orbital prob- 
lem, we analyze the point group symmetry of the two 
band model describing the dxz and dyz orbitals. We con- 
sider the behavior of the Hamiltonian under the tetrago- 
nal point group D^h = Civ '^ Ci where Ci is the inversion 
and C4u contains next to the identity E two four-fold ro- 
tations C4 one two-fold rotation C2, two mirror reflexions 
along the axis cr„ and two mirror reflexions along the di- 
agonals Gd- The Hamiltonian is invariant with respect to 



A rotation by 7r/2 that causes a sign change of an off- 
diagonal element of $ is therefore no indication for pair- 
ing in the d-wave channel. Thus assuming a rotation 
by 7r/2 (as generated by C4) yields a sign change of the 
off-diagonal element and no such change occurs for the 
diagonal element, we flnd C4$(fe^,fcj^) = ^{k^,ky), i-e. the 
gap belongs either to the irreducible representation Ai 
or A2. If furthermore the gap doesn't change sign upon 
reflection on the axis we conclude it is Ai, corresponding 
to s-wave pairing. This will be the case in our subsequent 
analysis of the numerical solution of spin fluctuation in- 
duced pairing. 
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Figure 2: (color online) Band occupation number along 
(0, 0) —f (tt, 0) —f (tt, tt) -^ (0, 0) for noninteracting case (red 
and blue) and with interactions (green and violet). 



C. Results 

In Fig. [2] we show the occupation number Up deter- 
mined from the full solution of the self consistent equa- 
tions in the normal state at T = 0.004eV. We com- 
pare our results with the corresponding occupation of 
the tight binding model without interaction at the same 
filling. The electron band closest to p = (7r/a, 0) under- 
goes a substantial distribution of carriers as it is being 
pushed very close to the Fermi energy. Similarly we ob- 
serve a decrease in the Fermi surface volume of the hole 
band centered around p = (0,0). Still the overall shape 
and topology of the various Fermi surface sheets are un- 
changed by many body interactions. 

In Fig. [3] we show the momentum dependence of 

the Qi — component of the effective interaction 
pai^a3,a4a2 (-p^ ^^ ^ Qy rj,j^-g -^ ^^^ ^j ^j^^ dominating 

components. Other matrix elements of Tph (q) have a 
similar momentum dependence. Finally Tph (q) and the 
particle particle interaction Tpp (q) behave very similar. 
The three panels show the effective interaction mediated 
by collective spin and charge fluctuations for three dif- 
ferent values of the Hund's coupling Jh- We clearly see 
that the effect of Jh is two-fold. On the one hand, larger 
values of the exchange coupling lead to an increase of 



the Stoner enhancement in Tpp and Tph- In addition, 
the effective interaction becomes increasingly more com- 
mensurate as Jh increases. The strong peaks close to 
p = (±7r/a, 0) and p = (0, ±7r/a) are consistent with the 
observed Bragg peaks for the magnetic ordering in the 
undoped parent compounds/^' 

In Fig. [4] we show the variation of the largest eigen- 
value A as function of the exchange and Hund's cou- 
pling J for two temperatures T = 0.004eV ~ 46K and 
T = 0.006eV ~ 70K. The enhancement of the effec- 
tive pairing interaction, discussed in Fig. 3, is the primary 
reason for the enhancement of the pairing strength and, 
in turn, of the leading eigenvalue A. We also flnd that 
A = 1 for J ~ 0.4eV, which would correspond to a crit- 
ical temperature Tc ~ 70K. While the above mentioned 
static approximation tends to overestimate T^ these re- 
sults demonstrate that experimentally relevant Tc values 
are clearly possible within the spin fluctuation approach. 
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Figure 3: (Color online) 0301,0204 = (0,0) component of 
the pairing interaction ^^311,0204 (p^,^^ = Q), Eq. iflSJI . for 
J = O.OeV (a), 0.25eV(b), and 0.50eV (c). Pairing interaction 
becomes increasingly commensurate as the Hund's coupling 
J increases. 
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Figure 4: (Color online) Largest eigenvalue A for the linearized 
version of Eq. dSI with T"^'^^'"'^"'^ (p,a;„ = 0) as a function 
of Hund's coupling J for T = 0.004eV and O.OOBeV. 



In Fig. [5] we show the momentum dependence of 
^xx (p) and /S.xy (p) as determined from the leading 
eigenvector of the Hnearized gap equation at T = 
0.006eV. The indicated diamond corresponds to the Bril- 
louin zone boundary, i.e. we plot the gap of the two xz 
orbitals within the unit cell in an extended zone scheme. 
The fact that both gap functions are of comparable mag- 
nitude reflects the fact that Cooper pairs are formed out 
of electrons in the same and in different d-states. The 
symmetry of the gap function is s-wave, i.e. it is in- 
variant with respect to the point group operations of the 
Hamiltonian. Simultaneous rotation of momenta p and 
orbitals by 7r/2 yields I^^^x {Px,Py) -^ \y {Py,-Px) and 
^xy{Px,Py) -^ -^yx{Py,-Px)- The latter expression 
explains the sign change of Aj;j, (p) upon rotation. It is 
a consequence of the s-wave symmetry in a two orbital 
problem where the xz and yz orbitals transform like the 
two dimensional coordinates. The fact that the diago- 
nal gap IS.XX (p) differs for momenta pointing along the 
two diagonals of the Brillouin zone is a consequence of 
the fact that the wave functions for the xz and yz or- 
bitals are different, see Ref.-li. Changing the value of 
the exchange interaction does not change the symmetry 
of the gap function. However, it signiflcantly affects the 
momentum dependence of Aaia2 (p)- As mentioned, the 
pairing interaction for small Jr is incommensurate with 
peaks rather far away from the ordering vector (tt/c, tt/o) 
of the antiferromagnetic state in undoped systems at am- 
bient pressure. On the other hand, for J = 0.25eV, the 
dynamic magnetic susceptibility and the pairing interac- 
tion Tpp (p) are peaked very close to (tt/o, tt/o). A com- 
mensurate pairing interaction can more efffciently change 
the sign of the gap function in momentum and orbital 
space, while incommensurations tend to frustrate an op- 
timally shaped pairing gap. This leads to the more com- 
plex pairing state for small J. 

Finally we determine the consequences of this gap func- 




Figure 5: (Color online) Momentum dependence of A^^ and 
Ajjy determined from the eigenvector corresponding to the 
leading eigenvalue of the linearized gap equation at T = 
0.006eV for Hund's coupling J = O.OSeV (top) and J = 0.25eV 
(bottom). Gaps of the two xz orbitals are shown in an ex- 
tended zone scheme. White diamonds indicate the Brillouin 
zone boundary. Red (light) and blue (dark) regions corre- 
spond to opposite signs of the gap. 



tion and analyze the gap anisotropy on the Fermi surface. 
From the self energy Sj^ {iijJn) we determine the quasi- 
particle energies £;;"'^ = E^P + S^'' (0) - ^l5a(3 and con- 
struct the quasiparticle energies of the superconducting 
state from the eigenvalues of 



Et Ar 



A, 



-E* 



(33) 



In Fig. [6] we plot the magnitude of the gap along the 
various sheets of the Fermi surface. The Fermi surface 
is constructed from the minima of the magnitude of the 
eigenvalues of hp. As shown in Fig. [6^, we flnd that 
in case of a small J the pairing interaction is more in- 
commensurate and the gap vanishes on line nodes on the 
Fermi surface. However, for larger J values we only flnd 
moderately anisotropic gap amplitudes on the Fermi sur- 
face, see Fig. [6b, c. The gap amplitude on the inner Fermi 
surface sheet around F is signiflcantly larger than the 
gap on the outer sheet, in agreement with recent ARPES 



experiments 



17.18 



This is a consequence of the fact that 



Axx (p) and Ayy (p) change sign close to the Brillouin 
zone center. The gap of the Fermi surface sheets cen- 
tered around M are considerably more anisotropic and 
could be responsible for the observation of anisotropic 
gaps.— 'i^>i^'i^>2£ In general, experiments that are sen- 
sitive to the minimum of the gap should therefore flnd 
much smaller typical gap values and more anisotropic 



a) 



<i 



b) 



< 






//i\ ft i( ;\-Ji. 
^1 






Figure 6: (Color online) Amplitudes of the gaps along four 
sheets of the Fermi surface for J — O.OeV (a), 0.25eV (b), 
and O.SOeV (c). While linearized gap equation cannot define 
the absolute amplitude Ao, the relative gap amplitudes are 
properly defined. 



gaps than measurements that are more sensitive to the 
largest gap values. 

Our calculation yields a fully gapped Fermi surface in 
the case where the pairing interaction is close to being 
commensurate. In this case the nodes of the gap are 
located between different Fermi surface sheets, explain- 
ing the dramatic change in the amplitude of the gap as 
one gets closer to the nodal lines (see Fig. 6). The po- 
sition of these nodes is not fixed by symmetry and, as 
is seen in case for more incommensurate pairing inter- 
actions, can in principle touch the Fermi surface (see 
Fig. 6 a). It is therefore an interesting question to ask 
what happens if one includes electron-electron overlap 
between different FeAs layers. This seems particularly 
relevant for the 122 materials where the outer sheet of 
the Fermi surface around F = (0, 0) increases its radius 
for increasing k^ M If the pairing interaction is predomi- 
nantly two dimensional, and determined by those Fermi 
surface sheets that are less dispersive in the z-direction, 
we expect that the position of the nodes is only weakly 
affected by the dispersion along k^. It is therefore eas- 
ily possible that at least one Fermi surface sheet touches 
the nodal plane for larger fcz values. The intersection be- 
tween nodal plane and Fermi surface would then yield a 
nodal line on the Fermi surface. This implies that one 
can easily explain fully gapped pairing states and states 
with line nodes with same pairing symmetry (s^) and due 
to the same pairing mechanism. Note, this is impossible 
for a d-wave pairing state, which will always yield line 
nodes given that the Fermi surface around the F point is 
closed. It is also impossible within a conventional s-wave 
pairing state where the sign of the gap is the same every- 
where. Thus, seemingly conflicting observations in dif- 
ferent FeAs-based systems do not necessarily imply that 
there are several distinct pairing mechanism at work. 

In summary, we determined the anisotropy of the spin 
fluctuation induced pairing gap on the Fermi surface of 
the FeAs based superconductors. For realistic param- 
eters we find a fully gapped state, while a measurable 
anisotropy remains for some Fermi surface sheets. This 
may explain the conflicting observations for the pres- 
ence of gap nodes obtained in NMR, penetration depth 
and ARPES experiments. It does explain the variation 
of the gap on distinct sheets of the Fermi surface, as 
seen in ARPES experimentsjii^ More generally, our re- 
sults demonstrate that a fully gapped superconducting 
state is fully consistent with an unconventional pairing 
mechanism. 
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